Two-soliton collisions in a near-integrable lattice system 
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We examine collisions between identical solitons in a weakly perturbed Ablowitz-Ladik (AL) 
model, augmented by either onsite cubic nonlinearity (which corresponds to the Salerno model, 
and may be realized as an array of strongly overlapping nonlinear optical waveguides), or a quintic 
perturbation, or both. Complex dependences of the outcomes of the collisions on the initial phase 
difference between the solitons and location of the collision point are observed. Large changes of 
amplitudes and velocities of the colliding solitons are generated by weak perturbations, showing that 
the elasticity of soliton collisions in the AL model is fragile (for instance, the Salerno's perturbation 
with the relative strength of 0.08 can give rise to a change of the solitons' amplitudes by a factor 
exceeding 2). Exact and approximate conservation laws in the perturbed system are examined, with 
a conclusion that the small perturbations very weakly affect the norm and energy conservation, but 
completely destroy the conservation of the lattice momentum, which is explained by the absence of 
the translational symmetry in generic nonintegrable lattice models. Data collected for a very large 
number of collisions correlate with this conclusion. Asymmetry of the collisions (which is explained 
by the dependence on the location of the central point of the collision relative to the lattice, and 
on the phase difference between the solitons) is investigated too, showing that the nonintegrability- 
induced effects grow almost linearly with the perturbation strength. Different perturbations (cubic 
and quintic ones) produce virtually identical collision-induced effects, which makes it possible to 
compensate them, thus finding a special perturbed system with almost elastic soliton collisions. 

PACS numbers: 05.45.-a, 05.45.Yv, 42.65.Tg, 47.53.+n 
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I. INTRODUCTION 

Soliton collisions is one of central topics of nonlinear- 
wave dynamics. In integrable systems, solitons are well 
known to emerge unscathed from collisions However, 
even small nonintegrable perturbations may render the 
phenomenology much richer, causing various inelastic ef- 
fects, such as trapping and formation of bound states, 
multiple-bounce interactions (where solitons separate af- 
ter multiple collisions) fractality in the outcome of 
the collision 0, 0|, and others. Such complex features 
are usually attributed to the excitation of soliton inter- 
nal modes |^,^|^, but more recently it was realized that 
they also occur due to the possibility for radiationless 
energy exchange (even in the absence of internal modes) 
between the colliding solitons, should the conservation 
laws allow it 6, 7J. The latter mechanism was both con- 
firmed by direct simulations of the corresponding nonin- 
tegrable models, and might be expected to ensue from 
the principle stating that any outcome compatible with 
the conservation laws may take place under appropriate 
initial conditions. 

Such effects suggest that the integrability is essentially 
tantamount to the strictly elastic character of the colli- 
sions , and warrant the importance of further studies of 
strongly inelastic collision effects produced by small con- 
servative perturbations added to basic integrable models. 
This general issue is of interest not only in its own right, 
but also for applications to nonlinear optical waveguides, 
as strong changes in the character of the interaction in- 



duced by a small perturbation may be naturally used 
in the context of switching, see e.g., @ and references 
therein. 

The objective of the present work is to consider such 
effects in collisions of dynamical (rather than topologi- 
cal) solitons in a discrete near-integrable system. As a 
matter of fact, the only integrable system which can be 
used in this case as the zeroth-order approximation is the 
Ablowitz-Ladik (AL) lattice @. Its well-known noninte- 
grable extension is the Salerno model (SM) 1^, which 
is produced by adding the integrability-breaking pertur- 
bation in the form of the onsite nonlinearity to the inte- 
grable AL system with the inter-site cubic nonlinearity. 
In order to test if the results that will be obtained below 
are generic, we will also consider an essentially different 
type of an integrability-breaking conservative perturba- 
tion, viz., the quintic onsite nonlinearity [its principal dif- 
ference from the cubic counterpart is that it breaks the 
integrability of both the AL lattice and of its continuum 
limit, i.e., the nonlinear Schrodinger (NLS) equation]. 

Thus, we introduce a general dynamical model based 
on the following equation: 



+ (l/2)(l-<5)|^„|2(V7„+i-f^„_i) 



e\^pn\'^^pn, 



Here ipn is the complex dynamical variable at the n-th 
site of the lattice, the overdot stands for the time deriva- 
tive, h is the lattice spacing, e is a real constant control- 
ling the quintic perturbation, and (5 is a real parameter 
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that accounts for the crossover between the AL {6 = 0, 
e = 0) and discrete-NLS {6 = 1, e = 0)^ Hmits. Equa- 
tion conserves two dynamical invariants, namely, the 
norm, 



Af=-^S2^n[l + h^{l-SM^\'] 



(2) 



and energy (Hamiltonian) , 



n 



E 



h^il-6)+e 



\n[^ + h^{^-S)\^|:n\''] 



■ 2(1 - 6) 



(3) 



While the discrete NLS equation, corresponding to 
6 = 1 and e = in Eq. (Q, has numerous physical re- 
alizations, the most important one being arrays of non- 
linear optical waveguides the AL model does not 
directly apply to many physical systems, because of the 
specific character of the nonlinear terms in it. However, 
a realization of the SM may be an array of strongly over- 
lapping nonlinear optical waveguides, especially the one 
following a zigzag pattern (similar to an array introduced 
in Ref. 13]). Indeed, the overlapping between adjacent 
cores will give rise, through the Kerr effect, to a nonlinear 
correction in the linear coupling between the cores, in the 
form of the terms ~ (1 — 5) in Eq. (^). It should be noted 
that, in this case, extra perturbation terms are expected 
too, such as (IV'ti-iP + IV'n+iP) V'ti (cross-phase mod- 
ulation). However, the results presented below clearly 
demonstrate that strong effects generated by small con- 
servative perturbations are essentially the same for dif- 
ferent perturbations, therefore we expect that taking into 
regard all the possible perturbation terms corresponding 
to the optical waveguides with strong overlap between 
the cores will not alter the results significantly. 

In some specific cases, soliton collisions in the SM have 
already been examined. In particular, a collision between 
a soliton and a reflecting wall, which is equivalent to a 
strictly symmetric collision between the soliton and its 
mirror image, were studied numerically in Ref. [T^ . One 
of our aims is to explore sensitivity of collisions to small 
asymmetries in initial phases and positions of the solitons 
in the actual two-soliton collision. Very recently, colli- 
sions in the (strongly nonintegrable) discrete NLS equa- 
tion were examined 15], and symmetry- breaking effects 
were found, along with sensitivity of the outcome to the 
location of the collision point. Here, we present results 
of collisions and their dependence on parameters in the 
model ^ with small 6 and e, i.e., close to the integrable 
AL limit. Together with the already available findings 
for the strongly nonintegrable case |l5| , they provide for 
a sufficiently comprehensive description of the coUisional 
dynamics of nontopological solitons in fundamental lat- 
tice systems. 



The AL model [Eq. (dJ with 5 = and e = 0] has 
exact soliton solutions of the form 



1 



sinh fi 



h cosh[fi{n — x(t))] 



exp {ik [n — x{t)] + ia{t)} 



(4) 

where the instantaneous coordinate and phase of the soli- 
ton are 



x{t) 
a{t) 



t sinh n 
xo + -TT (smfc) 



(cos k) cosh fi + k (sm k) 1 



5.) 



xq and ao are their initial values, while /i and k define 
the soliton's amplitude A and velocity V, 



A = h ^ sinh^ /i , V = (fih) "'^ (sinh /x) sin fc . 



(6) 



The infinitely long AL system has an infinite series of 
dynamical invariants, the lowest ones being the norm, 
lattice momentum, and energy. 



A^ = ^ln(l + /i2|^^ 



ih 



(V'nC 



■ ■0j^V'n+l), 
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(V'nV-'^^+l + V^V^n+l)- 



(7) 



(8) 



(9) 



Note that the norm of the general nonintegrable model 
(Q, given by the expression goes over into the norm 
in the limit 6 = e = 0, and the Hamiltonian (O 
of the nonintegrable model becomes, in the same limit, 
a linear combination of the norm Q and energy of 
the AL integrable system: H{6 = e = 0) = - {N + Q). 
It will be seen below that, as a matter of fact, the dif- 
ference between the exact norm and Hamiltonian of the 
full perturbed model with small 6 and e and those of 
the AL model is negligible. However, all the other dy- 
namical invariants of the AL model, including the lat- 
tice momentum 0, have no counterparts in the nonin- 
tegrable case. This is explained by the fact that each 
elementary dynamical invariant is generated by a certain 
continuum invariance of the underlying equation. In par- 
ticular, the norm and energy conservation are accounted 
for by the invariance against phase and time shifts, re- 
spectively, that remain valid in the nonintegrable system, 
while a hidden dynamical symmetry of the AL model 
which is responsible for the conservation of the lattice 
momentum is destroyed by the small perturbations. The 
momentum remains a dynamical invariant in continuum 
nonintegrable models (e.g., the NLS equation with the 
quintic term), but in the discrete setting it is conserved 
solely in the integrable case - obviously, a generic lattice 
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system is not invariant against arbitrary spatial transla- 
tions. Further consideration of this issue can be found in 
a recent work 16]. 

An issue related to the lack of conserved momentum 
is the (non)existence of exact traveling soliton solutions 
in nonintegrable lattice models. While the problem still 
waits for its full solution, several important theoretical 
results have been obtained. The inclusion of an appro- 
priate traveling ansatz in the discrete equation gives rise 
to a differential-delay equation whose steady states are 
the traveling wave solutions of the original differential- 
difference equation (see e.g., Refs. 0,0, as well as an 
earlier work by Feddersen, Ref. 0]). It is also pertinent 
to mention that moving solitons clearly persist in sim- 
ulations of perturbed systems, without any conspicuous 
loss, for long times, which is sufficient to study their col- 
lisions in numerical experiments without ambiguity (see, 
e.g., Refs. [11 and 13). 

If the given system differs from the AL model by small 
terms, a natural question is how strong actual destruc- 
tion of the former dynamical invariants, and especially of 
the momentum, which has a straightforward physical in- 
terpretation (and remains a virtually conserved quantity 
for free moving solitons, as explained above) will be in 
collisions between solitons. One of main objectives of the 
present work is to address this issue. We will conclude 
that the momentum conservation is strongly violated by 
the collisions, even if the perturbation parameters are 
quite small. 

For two broadly separated AL solitons with parameters 
fij and kj , the expressions Ql - © take values 



N. 



sol 



— 2 fij , Qsoi = 2 (sinh fij ) cos kj 



(10) 



Psoi = 2 (sinh fij) sin kj 

3 = i 



(11) 



The fact that only two exact and, plausibly, one approx- 
imate conserved quantities (the latter one is the momen- 
tum) constrain possible outcomes of the soliton-soliton 
collisions, which are characterized by two amplitudes, 
two velocities, and, in addition, may depend on the initial 
relative phase Aao = ao2 ~ chqi and positions (^0)2 and 
(xo)i of the solitons, suggests that the above-mentioned 
radiationless energy exchange between the two solitons is 
quite feasible. Furthermore, for slow solitons (small kj) 
the conservation of Q becomes an amplitude constraint 
[see Eqs. ()], and if amplitudes are small too (/ij 0), 
Q reduces to N, see Eqs. (jTUIl . so that there actually 
remains the single constraint in this limit case. 



II. NUMERICAL RESULTS 
A. Setting up the problem 

To perform simulations of the collisions, we notice that 
h can actually be scaled out from Eq. leaving S 

and e/h'^ as independent control parameters, therefore in 
what follows, we fix /i = 0.8. The parameters are varied 
in ranges corresponding to the weak quintic perturbation, 
e S [—0.01,0.01], and moderately weak Salerno's pertur- 
bation, 5 € [—0.08,0.08]. Equations (Q) were integrated 
by means of an implicit Crank-Nicholson scheme with 
the accuracy of 0((At) ) and with reflecting boundary 
conditions. The initial condition was taken as a super- 
position of two far separated solitons I^J that would be 
exact solutions of the AL model, and the numerical in- 
tegration was run until outgoing solitons were separated 
well enough. Their amplitudes and velocities after the 
collision, Ai, A2, Vi, and V2, were measured, and then 
the corresponding parameters jli, jl2, ki, and ^2 were 
found inverting Eqs. (P). We also checked to what ex- 
tent the dynamical invariants given by the unperturbed 
expressions (jlOII and Hll|) . as well as by the exact ones 
(PI and ((2Jl, were conserved. 

To present the results, we will focus on the symmetric 
collisions: fJ^i — IJ-2 — fJ', ki — ~k2 = k. While simula- 
tions were run for various values of the amplitudes and 
velocities, we display results for a case that turned out 
to be a typical one, adequately representing many oth- 
ers, for ^ = 0.75, k = 0.1. This implies Ai = A2 = 
A = 1.057, Vi ^ -V2 ^ V ^ 0.137. The initial phase 
difference Aao was controlled by setting agi = and 
choosing ao2 from the interval (— 7r,7r). The initial posi- 
tions of the solitons were taken as (a:o)i = —2^0 + Xc and 
(2^0)2 = 2^0 + Xc, with xo = 12; this provides for the large 
initial separation 2a:o = 24 between them, while Xc was 
chosen from the interval [0,1) to control the location of 
the collision point. 

The presentation of results is structured as follows: we 
first examine the effect of variation of the initial phase 
difference Aao and collision point Xc on the outcome of 
the collision. Then, we analyze how the approximate 
conservation of the expressions (|10() and Hll() correlates 
with the results. Finally, we examine the effect of com- 
bining the perturbation parameters e and 6, in order to 
demonstrate that the two perturbations may almost ex- 
actly cancel each other, thus making the collisions vir- 
tually elastic. We stress that results obtained at other 
values of parameters are completely tantamount to those 
displayed below, provided that e and i5 remain small. 



B. Sensitivity to the phase and position of the 
collision 

In Fig. n values of the soliton parameters after the 
collision are presented as functions of the initial phase 
difference Aao for the case of the SM perturbation. In 
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the part of the interval (— tt, tt) which is not included, the 
collision is almost completely elastic. It is obvious that, 
in the interval |Aq;o| ^0.5, the small perturbation is, in 
fact, a singular one, resulting in very strong effects (in 
other words, the elasticity of soliton collisions in the AL 
lattice is a very fragile feature). Note that, in the case of 
a weakly perturbed continuum NLS equation considered 
in Ref . ,6| , noteworthy inelastic effects in the collision of 
two solitons also took place at relatively small values of 
|Aao|. 



C. Dynamical invariants 

In the present case, the initial values of the expressions 
(PI and jni) are iV = 3, Q = 3.273, and P = (the 
absolute values of the initial momenta are |Pi| — IP2I = 
0.1642). Using the exact expressions Q and ©, we have 
checked that the net values of the norm and energy for 
the initial solitons and those observed after the collision 
are equal, in the case of the ordinary numerical accuracy 
employed, up to 10~* (in relative units); running simu- 
lations with higher accuracy (smaller At), it was possi- 
ble to check the norm and energy conservation with the 
accuracy of up to 10~^. Norm and energy loss due to 
radiation loss remained completely negligible in all the 
cases considered. As concerns the difference between the 
unperturbed expressions used in Eqs. (|10|l and the exact 
ones (0) and Fig. El demonstrates that the largest 
relative difference between them, which reflects a direct 
effect of the small perturbations, is AN/N ^ 10~^ for the 
norm, and AQ/Q ^ 3 x 10"'^ for the energy. However, 
the bottom panel in Fig. [21 shows that the momentum is 
not conserved in any approximation, in accordance with 
the fact that the perturbed system has no translational 
symmetry. 

The conservation of N suggests that a simple relation 
between the soliton amplitudes after the collision may be 
expected: according to Eqs. /ii -t- fj,2 must keep the 

original value with the accuracy ^ 10~^. On the con- 
trary, the momentum nonconservation promises a much 
worse accuracy in the prediction of a relation between 
the velocities. The conservation of Q does not provide 
for an essential additional information for small values of 
k (see above) , while large k implies the collision between 
fast solitons, when nontrivial effects will be very weak. 

The comparison between the actual results of the col- 
lision (dots) and predictions based on the approximate 
conservation laws for the values N and P of the two soli- 
tons in the form of Eqs. H10|) (dashed lines) is displayed 
in Fig. As is seen, the amplitude relation indeed fol- 
lows from the norm conservation in a very accurate form, 
while the conservation of the momentum may be traced 
in a very crude form only. 

Another salient feature of Fig. O is strong deviation 
of the dots from the diagonal point (1.054, 1.054) cor- 
responding to the values of the parameters before the 
collision. Such a feature was impossible in the case of 
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FIG. 1: Velocities and amplitudes of the solitons after the 
collision vs. the phase difference Aao, in the case S = 0.04, 
e = (the Salerno model without quintic terms). This and 
all the other cases are shown for fi = 0.75 and k = 0.1, see Eq. 
(HJ. Four different curves correspond to different positions of 
the collision point, Xc ~ 0, 0.25, 0.5, and 0.75. The collisions 
are strongly inelastic in the vicinity of Aao = 0. 



the collision of a soliton with its mirror image in the SM, 
examined in Ref. [T^ . A typical example of an inelastic 
collision (inducing this effect) is shown in Fig. 21 The 
major cause of the effect is the location of the collision 
central point Xc relative to the underlying lattice. 

Besides that, the phase difference between the collid- 
ing solitons may produce a similar symmetry-breaking 
effect. Indeed, if the AL solitons described by Eqs. Q 
and lO moving to the right and to the left [with k > 
and A: < 0, accordingly], are given phase shifts +A(f> and 
— A0, this is equivalent to the shift of the coordinate x, 
but solely in the expressions for the solitons' phases, by 
Ax — (t>o/k, which has equal signs for both solitons. This 
means the phase pattern of the two-soliton configuration 
gets shifted by Ax relative to the shapes of the collid- 
ing solitons, which is an obvious cause for the symmetry 
breaking. The fact that the above-mentioned position 
and phase factors do not affect the collision symmetry in 
the integrable AL model is another specific manifestation 
of its integrability. 
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FIG. 2: The collision-induced changes of the net norm, energy, 
and momentum for the two solitons, defined as per Eqs. ^ 
and vs. Aao for different values of Xc in the Salerno 

model (recall that the tilde refers to the post-collision values 
of the corresponding quantities). The quantities displayed 
in this figure are obtained by adding up their values for the 
two solitons (rather than by direct calculation for the whole 
system). If the norm and energy are defined by the exact 
expressions Q and rather than the approximate ones 
nun , they are completely conserved. 



D. The role of the perturbation strength 

In Fig. O one can observe that, for S = 0.04 and 
e = 0, the maximum possible sohton amplitude after 
the collision is Amax = 1.64, and the maximum possi- 
ble post-collision velocity is T^nax = 0.263. These val- 
ues [and, in particular, their deviation from initial ones, 
{A,V) = (1.054,0.137)] may be regarded as a measure 
of the departure of the perturbed model from the inte- 
grability. In Fig. [S] we use Amax and Knax to gauge the 
deviation from the integrable case with the increase of 
the perturbation strength (the cases of both the SM and 
quintic perturbations are shown). It is concluded that the 
weak perturbations generate quite large inelastic effects, 
the inelasticity increasing almost linearly with the per- 
turbation parameter. For instance, the value S — -1-0.08 
of the relative perturbation parameter in the SM model 
gives rise to collision-induced changes of the soliton's am- 
plitude by a factor of ~ 2, and of the absolute value of 
the velocity by a factor of ~ 3. 

Another noteworthy feature is the asymmetry of the 
plots in Fig. [S] The asymmetry is due to the fact that in- 
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FIG. 3: Relations between the solitons' amplitudes (a) and 
velocities (b) after the collision in the Salerno model. The 
dashed lines show the relations predicted by the norm, en- 
ergy, and momentum conservation for integrable AL chain 
Eqs. (ESJ and (|TTJ- Dots are numerical results for 2500 col- 
lisions, with values of the initial phase difference Aao taken 
from the interval [—1.25, 1.25] with a step of 0.01, and the 
collision-point's coordinate Xc taken from [0, 1) with a step of 
0.1. Note that the norm was taken in the approximate form of 
Eq. IIUL which pertains to the unperturbed Ablowitz-Ladik 
lattice. The relative nonconservation of this norm after the 
collision is ~ 10~^ (see the text), while the exact norm of the 
perturbed model, as given by Eq. ||5J, is conserved exactly, 
within the numerical accuracy. 



ternal modes in the colliding solitons can be excited only 
when S > (for e = 0) or when e < (for (5 = 0) [2l| . 
Hence, in these cases, we observe a combined effect of 
the radiationless energy exchange and the internal-mode 
excitation, while for i5 < 0, e = and 5 = 0, e > 0, only 
the former occurs. Naturally, the net nonintegrability- 
induced effects are stronger in the cases where the inter- 
nal mode can be excited. 
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FIG. 4: An example of a strongly inelastic collision between 
solitons in the AL system with a small on-site quintic per- 
turbation, 5 = and e = —0.01. In this case, Aao = and 
Xc = 0.2. 



(a) 



£=0 



-0.08 -0.04 0.00 0.04 0.08 

S 

2.5r 



(b) e=0 - 
/ 

/ 



-0.08 -0.04 O.OC 0.04 0.08 
S 




(d) S^O 



-0.010 -0.0C5 0.000 0.005 0.010 -0.010 -0.005 0.000 0.005 0.010 

e s 



FIG. 5: The maximum possible amplitude ylmax (filled cir- 
cles) and velocity Vmax (empty circles) of the soliton after the 
collision vs. the perturbation strengths 5 and e: (a,b) e = 
(the Salerno model); (c,d) (5 = (the quintic model). 



E. Compensation of perturbation effects 

In Ref. ||6|j, it was found that, in continuum models, 
inelastic effects in soliton collisions can be strongly sup- 
pressed if contributions from different perturbations can- 
cel each other. We have observed a similar feature in the 
present model. In particular, in Fig. |Slwe show the post- 
collision amplitude Ai versus Aao for three different per- 
turbations. It is clear that the cancellation takes place in 
the case (c), making the norm and momentum exchange 
an order of magnitude smaller than in the other cases. 



Similar compensation effects were observed for 5 = 0.08 
and £ — 0.01 in a wide range of soliton parameters, in- 
cluding non-symmetric collisions (between non-identical 
solitons). In fact, the possibility of the mutual compen- 
sation between the Salerno and quintic perturbation is a 
strong proof to the assertion that different conservative 
perturbations produce virtually identical effects, hence 
essentially the same results are expected from other per- 
turbations. 
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FIG. 6: The amplitude of the first soliton after the collision, 
Ai, vs. Aao for Xc = 0, 0.2, 0.4, 0.6, and 0.8, curves 1 to 5, 
respectively. The perturbations have (a) 5 = 0.08, e = 0; (b) 
S = 0,e = 0.01; (c) 5 = 0.08, e = 0.01. 



III. DISCUSSION AND CONCLUSIONS 

In this work, we have quantified properties of collisions 
between solitons in the Ablowitz-Ladik (AL) model with 
weak Hamiltonian perturbations. We have observed com- 
plex dependences of the outcomes of the collisions on the 
initial phase difference between the solitons and exact lo- 
cation of the collision point. Strong inelastic effects, in 
the form of radiationless energy and momentum exchange 
between colliding solitons, are generated by weak pertur- 
bations (for instance, a perturbation with the relative 
strength 6 = 0.08 gives rise to a change of the solitons' 
amplitudes by a factor exceeding 2). The effects pro- 
duced by different conservative perturbations are quite 
similar, suggesting that the results reported in this pa- 
per are generic. The exact and approximate conservation 
laws of the perturbed system were examined, with a con- 
clusion that the small perturbations very weakly affect 
the norm and energy conservation, but strongly destroy 
the conservation of the lattice momentum, which is ex- 
plained by the absence of the translational symmetry in 
nonintegrable lattice models. Statistical data collected 
for a very large number of collisions validate this conclu- 
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sion. Symmetry-breaking effects in the collisions (which 
are simply explained by the dependence of the result on 
the location of the central point of the collision relative 
to the lattice, and by the phase difference between the 
colliding solitons) were highlighted, and their magnitude 
was used to gauge the deviation of the perturbed model 
from integrability. It was also shown that, properly com- 
bining two different perturbations, it is possible to almost 
exactly cancel their integrability-destroying effects, thus 
constructing a perturbed system in which collisions arc 
practically elastic. 

In this paper, we were dealing with collisions between 
solitons with relatively large initial velocities. It would 
naturally be of interest to see how the picture presented 
is modified for smaller collision velocities, and, in partic- 



ular, to examine whether a fractal structure, similar to 
that observed in Ref. 0, can be found in the present 
model. This issue will be considered elsewhere. 



Acknowledgments 

We appreciate a discussion with M. Salerno. 
PGK gratefully acknowledges support from NSF-DMS- 
0204585. BAM appreciates hospitality of the Depart- 
ment of Mathematics at the University of Massachusetts 
(Amherst), and financial support from the European Of- 
fice of Aerospace Research and Development (US Air 
Force), provided through a Window-on-Sciencc grant. 



[1] M.J. Ablowitz and H. Segur, Solitons and the Inverse 
Scattenng Transform, SIAM (Philadelphia, 1981). 

[2] D. K. Campbell, J. F. Schonfeld and C. A. Wingate, 
Physica 9D, 1, (1983); M. Peyrard and D. K. CampbeU, 
Physica D 9, 33, (1983); D. K. Campbell and M. Peyrard, 
Physica D 18, 47 (1986). 

[3] P. Anninos, S. Oliveira and R. A. Matzner, Phys. Rev. 
D 44, 1147 (1991). 

[4] J. Yang and Y. Tan, Phys. Rev. Lett. 85, 3624 (2000). 

[5] P. G. Kevrekidis, Phys. Lett. A 285, 383 (2001). 

[6] S. V. Dmitriev, D. A. Semagin, A. A. Sukhorukov, and 
T. Shigenari, Phys. Rev. E 66, 046609 (2002). 

[7] S. V. Dmitriev and T. Shigenari, Chaos 12, 324 (2002). 

[8] A. B. Aceves, C. DeAngelis, G. G. Luther, A. M. 
Rubenchik, and S. K. Turitsyn, Physica D 87, 262 (1995); 
P. M. Ramos and C. R. Paiva, J. Opt. Soc. Am. B 17, 
1125 (2000); J. Meier, G. Stegeman, H. S. Eisenberg, Y. 
Silberberg, R. Morandotti, and J. S. Aitchinson, in Third 
IMACS meeting on Nonlinear Evolution Equations and 
Wave Phenomena, Book of Abstracts, p. 139. (University 
of Georgia: Athens, 2003). 

[9] M. J. Ablowitz and J. F. Ladik, J. Math. Phys. 16, 598 
(1975). 

[10] M. Salerno, Phys. Rev. A 46, 6856 (1992). 
[11] P. G. Kevrekidis, K. 0. Rasmussen, and A. R. Bishop, 
Int. J. Mod. Phys. B 15, 2833 (2001). 



[12] H.S. Eisenberg, Y. Silberberg, R. Morandotti, A.R. 
Boyd, and J.S. Aitchison, Phys. Rev. Lett. 81, 3383 
(1998). 

[13] N. K. Efremidis and D.N. Christodoulides, Phys. Rev. E 

65, 056607 (2002). 
[14] D. Cai, A.R. Bishop and N. Gr0nbech- Jensen, Phys. Rev. 

E 56, 7246 (1997). 
[15] I. E. Papacharalampous, P. G. Kevrekidis, B. A. Mal- 

omed, and D. J. Frantzeskakis, Soliton collisions in the 

discrete nonlinear Schrddinger equation, Phys. Rev. E, in 

press. 

[16] P.G. Kevrekidis, Physica D 183, 68 (2003). 

[17] G. looss and K. Kirchgassner, Commun. Math. Phys. 

211, 439 (2000). 
[18] G. Friesecke and R. L. Pego, Nonlinearity 15, 1343 

(2002). 

[19] H. Feddersen, in: Nonlinear Coherent Structures in 
Physics and Biology (ed. by M. Remoissenet and M. 
Peyrard, Springer- Verlag: Berlin, 1991), p. 159; see also 
D.B. Duncan, J.C. Eilbeck, H. Feddersen, and J.A.D. 
Wattis, Physica D 68, 1 (1993). 

[20] M. J. Ablowitz, Z. Musslimani, and G. Biondini, Phys. 
Rev. E 65, 026602 (2002). 

[21] Yu. S. Kivshar, D. E. Pelinovsky, T. Cretegny, and M. 
Peyrard, Phys. Rev. Lett. 80, 5032 (1998). 



